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ABSTRACT 



The correlation of meteorological events such as the jet stream, gravity waves and 
boundary layer circulation with the optical turbulence parameters, the transverse co- 
herence length r 0 and the isoplanatic angle 0 a is essential for interpreting and forecasting 
imaging and laser systems performance. In support of the United State Air Force Relay 
Mirror Experiment, the Naval Postgraduate School performed a series of six site char- 
acterization measurements near Kihei, Maui during August 1987 to July 1988. Spatial 
and temporal summaries of atmospheric events corresponding to the optical remote 
sensor data are presented using meteorological data from the National Weather Service 
Rawinsonde Observation Stations, synoptic charts, GOES-WEST infrared satellite im- 
ages and four Kihei, Maui rawinsonde datasets. To quantify the correlation between 
optical turbulence measurements and meteorological phenomena, four methods of cal- 
culating 0 T from rawinsonde data were investigated. Results show that existing 
rawinsonde systems arc inadequate for direct Cf calculation. However, moderate im- 
provements in the vertical resolution, the temperature resolution and probe response 
time will allow direct calculations of optical turbulence parameters from rawinsonde 
data. 
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I. INTRODUCTION 



Atmospheric turbulence produces optical phase and amplitude fluctuations that se- 
riously degrade laser beam propagation and the performance of imaging systems. In 
order to perform systems analyses of optical propagation, the spatial and temporal dis- 
tribution of atmospheric optical turbulence must be understood. Over the past two and 
a half years optical turbulence measurements in the form of isoplanatic angles, 0 0 , and 
transverse coherence lengths, r 0 , have been made at a variety of sites in the continental 
United States, as well as Maui. Hawaii and Vieques, Puerto Rico. Though a dominant 
application of these measurements involve laser technology and astronomy, the inter- 
pretation of the results are very valuable to the meteorologist. The optical data repres- 
ents an instantaneous summary of the optical turbulence along a vertical path. 
Imbedded within this vertical cross-section is information describing various synoptic 
and mesoscale phenomena such as the jet stream, gravity waves, low level jets and 
boundary layer circulation. By correlating meteorological measurements and optical 
turbulence parameters, the practical cause and effect relationships can be assimilated 
into prognostic tools. 

The United States Air f orce Relay Mirror Experiment (RME) requires a critical 
evaluation of atmospheric turbulence over the Maui. Hawaii experiment site. RME is 
designed to "assess how accurately an orbiting mirror can place (a) reflected (laser) beam 
on a ground target with low beam wander" (Nelson, 19SS). Observed beam jitter and 
wander arc a product of both the laser system and the atmospheric turbulence. Accurate 
measurements of the atmospheric contributions are vital for evaluating the optical sys- 
tem's performance. 

Between August 19S7 and July 19SS, the Naval Postgraduate School (NTS) per- 
formed six measurement sessions of approximately one week duration each in order to 
quantify the RME target board site's optical turbulence characteristics. NPS optical 
instruments were used to acquire atmospheric coherence length, r 0 , and isoplanatic angle, 
6 0 . measurements near sea level at Kihei, Maui. In addition to meteorological data from 
the National Weather Service Rawinsonde Observation Stations, synoptic charts and the 
GOES-WEST Infrared Satellite Images, four of the six Kihei data sets had on-site 
rawinsonde measurements. 
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This thesis attempts to correlate the RME meteorological measurements with the 
optical data. Initially, the National Weather Service measurements, charts and satellite 
images were used to compare meso- and synoptic scale events with the optical results. 
This method provided qualitative results. 

The major thrust of this thesis, however, materialized with the availability of high 
resolution rawinsonde soundings. These profiles finally provided the resolution neces- 
sary for attempting the direct computation of optical parameters from the rawinsonde 
data. The August and December 19S7 sessions used the VIZ WO-8000RP + rawinsonde 
system. With a vertical resolution of approximately 35 m, measurements were recorded 
up to 12 km. The April and July 19SS sessions made use of the Vaisala DigiCORA MW 
1 1 system. This system profiled meteorological data up to approximately 24 km with a 
sampling rate of 2 seconds (surface to 2 km) and 10 seconds (above 2 km). The S m and 
40 m vertical resolutions of the DigiCORA system were considered reasonable for com- 
puting the optical turbulence directly from the rawinsonde data. 

The optical variables measured. /•„ and 0 0 , depend on an integral of the index of re- 
fraction structure parameter. Q. along the vertical path. The parameter Q is propor- 
tional to the temperature structure parameter O t which can be computed from the 
rawinsonde data, four methods for computing the temperature structure parameter Q 
from the rawinsonde data were investigated: 

1. the Simple Potential Temperature Differences (SPTD) liquation, 

2. the Potential Temperature Differences with a Running Mean. 

3. the Bulk Gradient Equation (BG) and 

4. conditional data sampling based on the Richardson Number. 

Each Q method was tested by calculating r 0 and 6 0 from the Q profile. These results 
were ultimately compared to the optical measurements simultaneously acquired over the 
same site. 



II. BACKGROUND 



A. HISTORICAL REVIEW OF CORRELATION RESEARCH 

International scientific interest in the meteorological influence on "seeing", a rating 
system based on the angular size of a star image, flourished in the mid to late 50's. 
Following a brief decline, interest was reestablished in 1962. The new upper atmospheric 
and outer space environmental data introduced by artificial satellites and space probes 
provided a catalyst for the development of atmospheric physics. One of the dominant 
questions addressed by the melange of scientific backgrounds required in the atmo- 
spheric physics field involved how meteorological conditions aflcct "seeing" (Kucherov. 
1966a). 

In October 1962, the issue of observatory site selection and related astroclimatic 
topics were discussed at the International Astronomical Union Symposium held in 
Rome. Geophysicists, meteorologists and astronomers worked in cooperation to study 
the interrelationships between atmospheric conditions and astronomical "seeing". Par- 
ticipants resolved that placing telescopes at certain heights above the ground would im- 
prove image quality (Kucherov, 1966a). 'lhe following summarizes some of the 
pertinent articles presented at the Third All-Union Conference held four months earlier 
in Kiev. 

Kucherov (1966b) observed "shimmer angle" and meteorological conditions at one 
seashore and two inland sites over a period of approximately two years (1960-1962). 
The results of that study confirmed the relationship between synoptic fronts (mostly cold 
fronts) and astronomical "seeing". It should be noted that 25% of the cases display no 
relationship. 

Bol'shakova and Darchiya (1966) investigated the causes of "seeing" condition 
changes. While they acknowledge the synoptic fronts as an influence on "seeing", they 
suggest that "seeing" is affected by frontal passages only when air masses with optical 
characteristics differing from the original advect over the site. Using statistics, they 
support their claim. To a first approximation, 50% of the cases show the synoptic front 
causing a deterioration in "seeing", or alternatively, absence of a synoptic front ensures 
high "seeing" conditions. The rest of the cases display the reverse situation. In their 
opinion, the analysis must be repeated taking into account the type of front, the indi- 
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nation of the frontal plane, the direction of the front's motion, and ideally, other detailed 
meteorological information defined in various atmospheric layers. 

Zinchenko (1966) comes to the following conclusions: 

1. The lower-atmospheric turbulence is of decisive significance in all that concerns 
"seeing": 

a. "optical inhomogcneitics are most effective in the ground layer, where the air 
density is highest:" 

b. "turbulence intensity is greatest in this layer"; and 

c. "eddies of minimum size occur in the ground layer". 

2. Contradictor} - results exist with the applicability of Richardson number(Ri) as a 
criterion of high atmospheric turbulence. 

The main reason supporting the latter point is that turbulence intensity is not a 
single-valued function of Ri, though a dependency on Ri exists. Zinchenko explains, 
when Ri < 1 , one can qualitatively accept that turbulence set in: and, the amount of ac- 
cumulated eddy energy in steady-state conditions increases with the decrease in Ri. 
(Note: Current understanding requires Ri<0.25 for the above to be true.) Though lack 
of data prohibited an exact definition of the layer directly responsible for the blurring 
of the stellar diffraction disks. Zinchenko found that turbulence docs show a tendency 
to be restricted to the 0. 5-2.0 km layer. 

B. THE RELAY MIRROR EXPERIMENT 

The Relay Mirror Experiment (RME) is a United States Air Force Experiment 
monitored by the Kirtland Air Torcc Weapons Laboratory. Albuquerque, New Mexico. 
RME is designed to "assess the ability to place a laser beam on a ground target with low 
beam wander" (Nelson. 19SS). As described by Nelson. RME involves the propagation 
of three ground to space (uplink) laser beams as well as the reflection from space to 
ground (downlink) of one of these beams (see figure 1). Two of the uplink lasers, one 
from AMOS (Air Force Maui Observatory Station), a site located at the Haleakala 
crater lint, and the other from Kihci. a site located at the volcano's base, serve to align 
a relay mirror located on a payload orbiting over the Maui. Hawaii stations. With the 
mirror aligned, a 1.06 n m laser travels from AMOS to payload and is relayed down to 
a hexagonal target board at Kihci. The telescopes and optical detectors on the target 
board record the centroid motion and beam size of the reflected 1.06 ,u m laser beam. 
Atmospheric turbulence has an important effect on both uplink and downlink beams 
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(Nelson. 1988). Understanding and characterizing the atmospheric turbulence over the 
Hawaiian experiment site is crucial for assessing the success of the RME project. 

C. HAWAIIAN CLIMATOLOGY 

The following is a brief summary of Hawaiian climatology with a focus on the main 
experiment site in Kihei, Maui. 

The tropical island of Maui is located at 156° 15' longitude and 20° 45' latitude 
(Blumenstock and Price, 1967). The island is the second most southern island in the 
Hawaiian chain and has an extreme length of 77 km, maximum width perpendicular to 
extreme length of 42 km and a maximum elevation of 3.1 km. Nicknamed "The Valiev 
Isle", this 1S86 square km island is comprised of an isthmus linking two volcanos. The 
older volcano, Pui Kukui, is in the northwest sector reaching a peak of 1.8 km. To the 
southeast is the dormant volcano. Haleakala. Though Haleakala defines the island's 
highest point, it's crater is one of the largest on earth: 12 km long. 4 km wide and about 
900 m deep (Creal et al., 19S7). Topographically induced high winds crossing the sum- 
mit are common. 

1. Climatic Regions of Hawaii 

The tropical climate on Maui varies greatly with location. The dominate sub- 
climates include the following: 

1. Rainy MOUNTAIN slopes on the north, northeast WINDWARD SIDE (above 
600 m); 

2. Moderately rainy and mild WINDWARD LOWLANDS (below 600 m) to the 
north, northeast; 

3. Dry LEEWARD LOWLANDS to the south, southeast: 

4. Sometimes rainy, sometimes dry INTERIOR LOWLANDS between the two 
volcanic peaks; 

5. Dry. cold HIGH MOUNTAINS above 600-900 m (Blumenstock and Price, 1967). 

Rainfall measurements reveal the sharp weather contrasts over this island. On 
the leeward coastal areas such as Kihei, and mountain summits, average annual rainfall 
is 0.5 m. The lower windward slopes of a high mountain and at near the summit of 
lower mountains (western Maui) have an average annual rainfall of 7.6 m (Blumenstock 
and Price. 1967). 

2. Circulation 

The dominant atmospheric circulation for the llawiiain islands is the north- 
northeasterly trade winds. The prevalence of this flow reflects the seasonal north-south 
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movement of the North Pacific High. From May to September (summer), the North 
Pacific High is to the north. During this season, the southern sectors of the counter- 
clockwise flow around the High is such that trades prevail SO-95% of the time. From 
October to April (winter) the North Pacific High drops south such that the heart of the 
trades pass well to the south of the islands. The frequency of trades during this season 
drops to 50-S0% (Blumenstock and Price, 1967). 

Land and sea breeze circulation occurs on a much smaller scale. Kihei is unique 
in that land and sea breezes are present even during moderate to strong trades 
(Blumenstock and Price, 1967). The sometimes brisk sea breezes occur from forenoon 
to early evening. During the night, a very gentle land breeze flows until shortly after 
sunrise. The Kihei region also receives a nighttime Katabatic (downslope) flow compli- 
menting the land breeze already in progress (personal observation). 

3. Temperature, Humidity, Cloud Cover 

Hawaii's mean annual temperature range generally fluctuates less that 6 degrees 
Celsius (Blumenstock and Price. 1967). Daily summer temperatures at 2.5 km are be- 
tween -1C (nighttime) and IOC (daytime). The winter range at this height is 0-10C. In 
the lowlands, daytime temperature is between 21-27C, while nighttime is in the 16-2 1C 
range. The warmest months in the dry areas such as Kihei are August and September 
(Blumenstock and Price. 1967). 

Generally humid conditions exist with high clouds except in the driest coastal 
regions and at high elevations. In the absence of synoptic scale events. Kihei begins a 
day with locally clear skys. A backdrop of clouds clinging to the north and west slopes 
of I laleakala is common. Throughout the day Kihei experiences a growing and decaying 
of cloud cover. From sunset to approximately 2000 hours local time, the partly cloudy 
sky cover over the site dissipates. Because of the dramatic climate changes over such a 
small island, clouds are generally visible at all hours from the Kihei experiment site in 
some section of the sky (Blumenstock and Price, 1967). 

4. Storms 

Generally. Hawaiian weather falls into trade wind and non-trade wind condi- 
tions. During trade wind conditions one can expect trade wind showers on the E and 
NE (windward) coasts; possible light trade wind showers over the S and S\V (leeward) 
areas; and, sea breezes ascending mountain slopes potentially producing light to occa- 
sionally heavy showers in the late afternoon and evening. Non-trade wind conditions 
are associated with light and variable winds; moist air slowly advecting onto the islands 
from the southeast, south or southwest; intense local storms; and, the possibility of 



major and minor storm systems passing through the islands (Blumenstock and Price, 
1967). 

Major storms for Hawaii fall into one of four classes of disturbances: 

1. Cold Front Storms - winter storm systems associated with cold fronts that ap- 
proach from the northwest end of Hawaii, and move eastward barely reaching 
Oahu. 

2. Kona Storms - winter storm systems similar to cold front storms; winds are gener- 
ally from "kona" or leeward directions relative to the trades; well-developed kona 
storms are more widespread and more prolonged than the usual cold front storms. 

3. Tropical Storms - similar to hurricanes but with winds below 120 km hr. 

4. Hurricane - the least prevalent storm system: in 63 years only four hurricanes have 
crossed over the islands. Unlike the Cold Front and Kona Storms, the Tropical 
Storm and Hurricane can occur at any time. They are most likely to occur, how- 
ever. between the months of July to December (Blumenstock and Price. 1967). 

These events will be discussed in greater detail as they pertain to a specific case study. 



S 



III. ATMOSPHERIC OPTICS THEORY 



A. ATMOSPHERIC INFLUENCE ON LIGHT PROPAGATION 

Light propagates through a medium in the form of a wavefront, "a surface over 
which an optical disturbance has a constant phase" (Hecht and Zajac, 1974). Fermat's 
Principle explains that the route light travels will be the smallest Optical Path Length 
(OPL). This is expressed mathematically as (Hecht and Zajac, 1974) 

OPL = [' n{s)ds, (1) 

where n is the index of refraction (the "speed of an electromagnetic wave in vacuum to 
that in matter"), and s is the path length. The phase of a propagating wave will be 
k = 2-/7 times OPL (Eq. (1)). Optical turbulence arises when the accumulated phase 
along adjacent paths is different. Phase distortion in the atmosphere occurs when a 
wavefront encounters random irregularities in the atmosphere's index of refraction. 
Such irregularities are due to the turbulent mixing processes of the atmosphere. Thus, 
when a wavefront from a laser propagates vertically through the atmosphere, accumu- 
lated. random phase differences will degrade the laser beam and image system perform- 
ance. Typical atmospheric distortions result in: laser beam centroid wander, 

scintillation, image breakup and blurring. 

The study of these atmospheric random turbulence processes requires the utilization 
of statistical quantities, means, variances and high moments, f requently, variances and 
the higher order moments are not single-valued, unique or even bounded. Also, the at- 
mosphere is neither homogeneous nor isotropic. Kolomogorov and Tatarski introduced 
the concept of structure functions and local homogeneity to resolve these problems 
(Tatarski, 1961). 

Kolmogorov and Tatarski define the structure function. D R {T) , for an arbitrary 
spatially distributed random variable R(r) as 

D R (r) = <{R(F ] )-R(7 2 )} 2 >, (2) 

where F = F, — F 2 , and the angled brackets indicate an ensemble average of the enclosed 
quantity. Assuming the random process generating R is homogeneous and isotropic, 
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then only the magnitude of r is important, r = |r, — f 2 |. Local homogeneity assumes that 
the mean </?(/•,)> = </?(r 2 )> (Clifford, 197S). Restricting the separation r to the inertial 
subrange of turbulence, Kolmogorov (1961) develops a universal form of the single 
structure function 



where Q is a structure parameter, a measure of "the total amount of energy in the tur- 
bulence" (Clifford. 1978); r is the distance between the two sampled points; and, r must 
be between the turbulence innerscale, / 0 . and the turbulence outerscale, L 0 . Eddy sizes 
below the innerscale correspond to the eddy's energy dissipation through viscous effects. 
The outerscale corresponds to "the largest scale size for which the eddies may be con- 
sidered to be isotropic" (Clifford, 1978). 

For optical propagation through turbulence, the characterizing quantity is the index 
of refraction structure parameter (Walters and Kunkel, 1981) 



where = n{i\) and n 2 = n(r 2 ) arc the atmospheric index of refraction at points and r 2 
respectively. Though Q is difficult to measure directly, profiles calculated from 
rawinsonde sensor measurements show that Q generally decreases with altitude (Nelson, 
19SS). A more practical approach for measuring atmospheric turbulence is to measure 
the temperature structure parameter Q. 



Since the atmospheric index of refraction 



depends on pressure I* and temperature T (Tatarski. 1961), Q is related to Q by sub- 
stituting Eq. (6) in Eq. (4) to get 




( 3 ) 



n- 1 = (79.vl0 ~ 6 ) -jr 



( 6 ) 



Cl = (79 x l(T 6 -p-) 2 Cf. 



( 7 ) 
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assuming turbulence is isobaric. Specific methods for calculating Q will be discussed in 
Chapter 5. 

B. OPTICAL VARIABLES 

As described earlier, the laser beam is affected by atmospheric turbulence as a result 
of index of refraction fluctuations along the turbulent propagation path. Two parame- 
ters that characterize the optical turbulence generated by the atmosphere are: the 
transverse coherence length, r 0 , and the isoplanatic angle, 0 0 . 

1. Transverse Coherence Length 

The transverse coherence length, r 0 , is a measure of the lateral autocorrelation 
length of the electric field propagating through the turbulent atmosphere. Because of 
the Fourier transform properties of an imaging system, r 0 is also a measure of the spatial 
frequency response of the atmosphere, the atmospheric modulation transfer function 
(MTF). Viewed as a component of a linear system, the atmosphere acts like a low pass 
filter suppressing the high spatial frequencies, or details of an image. 

'file perturbation theory developed by Tatarski and extended by Fried (1966) 
and Lutomirski and Yura (1971) shows that the long term, modulation transfer function 
of the atmosphere is expressed in terms of the wave structure function D by 

MTF(v) = ex P ( j-!-). (8) 

where v is the spatial frequency in cycles m. / is the wavelength of light and /is the focal 
length of the optical system. For a plane wave, L) has the form 



l*L 



D{p) = 2.91 A’V /3 



Q(z>/r. 



*•0 



(9) 



where k is 2rr//. the wavenumber; p = |F, — r 2 | = 7/v, a distance measured at the entrance 
pupil of the optical system. L is the optical path length and Cj; is the index of refraction 
structure parameter. 

The general format of MTF(v) is exp( — ap ! 3 ) where a is a turbulence and wave- 
length scaling parameter. MTF(v) is frequently written as: 

.1/7TI, ■) . exp | - 3.44( f>\ - c.v,,| _ 3.44( / ) s ' 3 ], 
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representing Fried's (1966) version of the long term atmospheric modulation transfer 
function. Equivalently. Lutomirski and Yura's (1971) mutual coherence function (MCF) 
is 



M CF(p) = exp [ — ( —■ ) 5 ^ 3 ], (11) 

where p 0 is an autocorrelation length of the electric field. 

The parameters r 0 and p 0 are different measures of the electromagnetic (EM) 
wave's coherence as restricted by the atmospheric turbulence. They are related by 

r 0 = 3.44 3/5 p 0 = 2.1p 0 . (12) 

Using Eq. (S) and (9). the coherence length for a plane wave, r 0 becomes 



pl 



r 0 = 2. l[1.46/< 2 



Cn{2)dl-] 



■--,-3/5 



(13) 



The parameter r g represents the distance traverse to the direction of propagation where, 
on average, the correlation of the electric field is c -3 " (f ried. 1966). 

To interpret r 0 consider the following: 

1. For an optical aperture d smaller than the coherence length, the electric field across 
the aperture is coherent. The angular resolution will be proportional to /./«/ . 

2. For an optical aperture larger than the coherence length, the correlation across the 
aperture is limited to regions of r 0 or p 0 (depending on how the correlation is de- 
fined). 'l he EM intensities across the aperture will add (energy is conserved). The 
optical system's average angular resolution will be proportional to /.jr 0 (Walters, 
Favier and Hines. 19 7 9). 

3. Visually, for apertures smaller than /•„. atmospheric turbulence results in image 
centroid motion, a process called, "tilt". In telescopes larger than r 0 , the image is 
broken up into multiple spcckels, each associated with a r 0 sized region. The angle 
z./r 0 determines the overall envelope of the complex image. 

4. With respect to measuring the atmospheric optical turbulence, large r 0 values, such 
as 2O0-4UO mm, indicate a small amount of optical turbulence is present along the 
integrated path. Conversely, small r 0 magnitudes, such as 20-40 nun. imply a large 
amount of optical turbulence exists along the integrated path. 

2. Isoplanatie Angle 

The isoplanatie angle, 0 0 , is an angular measurement of spatial coherence be- 
tween two intersecting rays (Walters. 1 9 S 5 ) . Given two light sources a small angle 0 



apart, the random atmospheric phase fluctuations along the two paths are different. 
Figure 2 displays rays /•,, and >- 31 interacting with three different turbulence patches. At 
the atmosphere top the two paths undergo unique path distortions. At the surface (near 
point 1). the turbulent interaction is common to both rays. 6 0 is the angle between rays 
/•,, and r 31 such that the accumulated phase fluctuations are correlated to within c _1 of a 
perfect correlation. 

Fried (19S2) develops a mathematical representation for 6 0 for apertures d > > 

>o : 



0 0 = [2.91/S 



Q(z)z^dz}- 



^0 



(14) 



where k is 2 r//, wavenumber of light: C- is the refractive index structure parameter: and. 
z is the altitude from the giound. 

Typical Q profiles display a rapid decrease with altitude. The r 5 3 term is sig- 
nificant in that it emphasizes the optical turbulence at higher altitudes. The net product 
of these two terms (the integrand) results in a profile of relatively constant values from 
a few kilometers to 10— 15 kilometers (Walters. 1985). 

The practical interpretation of 0 0 is similar to r 0 . When 9 0 magnitudes arc small 
(1-2 a rad), the optica! turbulence along the integrated path and especially around the 
tropopause level is large. Conversely, when values are large (14-20 ,u rad), optical 
turbulence is small. 




Fig. 2. Diagram of lsoplanatic Angle, 0 0 . 
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IV. EXPERIMENT/DATA 



A. RME KIHEI DATA COLLECTION 

The initial sampling of isoplanatic angles, 8 0 , took place on the summit of \1t. 
Haleakala in Maui, Hawaii. These measurements commenced on 17 September 19S5 
and concluded on 1 February 19S6. Simultaneous meteorological surface measurements 
at the site provided some understanding of localized atmospheric events. A year and a 
half later, both 6 0 and r 0 were measured at the base of Haleakala (at approximately sea 
level), in Kihei. Six sampling sessions of four and five nights took place between August 
19S7 and July 19SS. Two Hawaii Rawinsonde Observation stations, Hilo and Lihuc, 
provided raw data listings for each of the 19S7-198S cases. Additional on-site 
rawinsonde launches supplemented the optical data for four out of the six sessions (Ta- 
ble 1 ). The meteorological data for the August and December 19S7 experiments used 
the VIZ WO-SOOORP + rawinsonde system with 10-sccond vertical resolution. The 
April and July 19SS experiments employed the Yaisala DigiCORA MW 11 system with 
2- and 10-sccond vertical resolution. All four cases used the OMEGA Navigational 
Network for wind measurements. Intermittent GOES- WEST Infrared Satellite Images 
complimented the rawinsonde data for all the cases. The primary focus of this study 
will be on the April and July 19SS experiments, due to the higher vertical resolution 
provided by the DigiCORA dataset. 



Table 1. MAIN OPTICAL AND METEOROLOGICAL RME SENSORS 



EXPERIMENT DATES 


OPTICAL DATA 


RAWINSONDE SITES 




0 o 

(a rad) 


I'o 

(min) 


Hilo. 

I Iawaii 


Lihue. 

Kauai 


Kihei. 

Maui 


17 Sep S5-1 Feb 86 


X 










1S-22 Aug 87 


X 


X 


X 


X 


X 


24-27 Oct S7 


X 


X 


X 


X 




5-9 Dec 87 


X 


X 


X 


X 


X 


27 Feb- 1 Mar 8S 


X 


X 


X 


X 




2~ Apr- 1 May 88 


X 


X 


X 


X 


X 


7-11 July SS 


X 


X 


X 


X 


X 
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B. KIHEI RME EQUIPMENT 

1. Optical Equipment 

Optical measurements were taken by an isoplanometer and transverse coherence 
length sensor designed built by Dr. D.L. Walters, NTS. Tor specific details on each in- 
strument the reader should pursue Stevens (1985) for the isoplanometer and Walters, 
Favier, and Hines (1979) for the transverse coherence length sensor. Isoplanatic angles 
were measured by an isoplanometer mounted on a portable 200 mm Celestron telescope. 
The software was run on an HP 217 computer. Isoplanometer data were automatically 
sampled once every second and displayed real-time on a time verses isoplanatic angle 
plot. Before storing the data, the samples were converted into 10 second averages. The 
complete isoplanatic angle system is capable of operating independent of human inter- 
action for 1-3 hours provided polar alignment is perfect and the chosen star's zenith 
angle does not exceed 40 degrees. 

The transverse coherence length sensor was mounted onto a portable 355 mm 
Celestron telescope. Software used for real-time processing required an operator to 
evaluate the image and stellar positioning before each sample. This added interface with 
a trained researcher served to suppress instrumental artifacts. The average sampling rate 
was between 1-2 minutes sample. Processed data were stored on an IIP 300 series 
computer. 

2. Meteorological Equipment 

Meteorological measurements specific to this thesis were a product of the 
Yaisala DigiCORA MW 11. This solid state, digital design rawinsonde package inte- 
grates the receiver, data processors and operator display into one unit. High speed 
precision sensors permit rapid sampling (at approximately 2 second intervals) resulting 
in high density data acquisition. The 95 mm by 145 nun (approximate size) RS 80 
Rawinsonde (Navaid Radiosonde) package was launched on an unravelling tether at- 
tached to a 200 gm balloon filled with helium. Pressure values were measured by a 
BAROCAP pressure sensor, a steel capacitive aneroid with a resolution of 0.1 hPa and 
an accuracy of 0.5 hPa. Temperature sampling was taken by a THERMOCAP sensor, 
a small capacitive bead in a glass casing. Resolution of the T1IERMOCAP was 0.1C: 
accuracy was 0.2C. The HL'MICAP, a thin film humidity sensor, measured humidity 
with a resolution of 1% relative humidity (RII) and an accuracy of 2% R1I. According 
to the Yaisala Ref. A0435. this sensor produced a reliable response even in low temper- 
atures and after exposure to condensation. The measuring range published was 0-100% 
Rll (Yaisala. 1986a). 
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The commutation between the sensors and the transducer electronics was ac- 
complished by a digital electronic switch. A water activated battery provided the RS 80 
"sonde" with a light weight power supply capable of keeping the package activated for 
the entire surface to 30 km (balloon permitting) launch sequence (Vaisala, 1986b). 

The telemetry frequency was between 400-406 MHz. Wind data for this exper- 
iment were calculated from the OMEGA Navigational System. The RS 80 received 
signals from the OMEGA Navaid Network (11-14 kHz), superimposed this information 
onto the thermodynamic data and relayed the sampling to the DigiCORA MW 11 
(Vaisala, 1986b). 

C. KIHE1 SONDE DATA PREPARATION FOR ANALYSIS 

The initial form of the DigiCORA MW 11 rawinsonde output used in this research 
began as a hardcopy listing of time in minutes and seconds, height in meters, pressure 
in millibars, temperature in Celsius, relative humidity in percent, dewpoint temperature 
in Celsius, mixing ratio in grains per kilogram, wind speed in knots and wind direction 
in degrees. Using a Dest PCSCAN scanner, the powerful IBM mainframe XEDIT. a 
software package capable of transfering IBM files to a IBM PC-AT, and the instant 
graphics display of QL'ATTRO. the initial hardcopy was duplicated as an ASCII file. 

Thermodynamic sampling rate for the first two kilometers was every two seconds. 
Beyond this threshold, sampling continued once every 10 seconds. Due to the limita- 
tions in the Omega Navigational system, winds were measured even - 10 seconds 
throughout the launch. A linear interpolation between the 1 0-second wind measure- 
ments was applied to the initial two kilometers remedying the inconsistent sampling rate. 

Variables calculated from the raw data include potential temperature, potential 
temperature gradient, shear. Richardson number. O r . Q. r 0 and 6 0 . 
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V. ANALYSIS/RESULTS 

To strengthen the link between the measured optical parameters and the 
meteorological data collected, a number of methods for calculating r 0 and 0 0 values were 
pursued. As explained earlier, both r 0 and 6 0 are a function of Q , which, in turn, is a 
function of O t (see Chapter 3). In this chapter, four methods for calculating O t from 
rawinsonde data will be investigated. However, the first step is to define an appropriate 
conservative temperature quantity to be used in the O t calculation. 

A. STANDARDIZING THE TEMPERATURE PROFILE 

An air parcel vertically displaced in the atmosphere experiences a temperature 
change as the parcel's pressure equalizes with the ambient environment (Clifford, 1978). 
Permitting the air to be moist but unsaturated, and assuming the atmosphere is in 
hydrostatic equilibrium, the rate at which the lifted parcel's temperature will decrease is 
the dry-adiabatic lapse rate. This is approximately 9.767C km. If the parcel is saturated, 
the release of latent heat produces a lapse rate that is less than the dry-adiabatic lapse 
rate (Huschkc. 1970). For this study, it is important to work with conservative variables, 
a variable that is independent of position in space. As described above, temperature 
does not fit this description. The conservative quantity (independent of position) best 
suited is potential temperature (Clifford. 197S). By definition, potential temperature is 
"the temperature which the parcel of air would hat e if it were expanded or compressed 
adiabatically from its existing pressure and temperature to a standard pressure />„ (gen- 
erally taken as lOOOmbf (Wallace and Hobbs. 1977). In short, the temperatures sam- 
pled vertically are standardized to one pressure level. Mathematically, the potential 
temperature expression is 



Pa -B- 

0 =T(-y)c r , (15) 

where 0 is potential temperature in kelvin. T is temperature in kelvin. 1\ is a standard 
pressure such as lOOOmb. P is pressure in millibars. R is the gas constant for dry air and 
C f is the specific heat of dry air at constant pressure (1.005 J g K) (Huschke. 1970). 
Note: Following the usual convention. Q shall hereafter be written in place of Q 

(Clifford, 19~$). 
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B. CALCULATING C> : FOUR CASE STUDIES 

1. Case I: Simple Potential Temperature Differences 

Case I takes the given mathematical definition, 

.2 <[(02-0 2 )-(©l-0,)] 2 > 

C 7 -= — (16) 

r 

and assumes © 2 = 0, . Physically, this presumes homogeneity between neighboring 
samples. In other words, the samples arc taken within the inner and outerscale of the 
turbulent layer (Tatarski. 1961). Note: the angled brackets indicate an ensemble average 
of the enclosed quantity. Associated with the above and all the cases hereafter is the 
assumption that the sampled layer is isotropic (Tatarski, 1961). 

Unlike most Q- calculations, the dataset employed by this research consists of 
vertical temperature samples. For the above assumption to be true, the rate of sampling 
must be sufficiently high with respect to the potential temperature gradients. If the 
sampling rate is too slow, trends in potential temperature will show up in the differences, 
(0 : — 0,). Figure 3 displays such trends first between 2 and 5 km then again above 1 2 
km. Note: The large increase in potential temperature differences in the stratosphere 
was found in all the profiles examined. 
a. Result!, 

Using the entire potential temperature profile for the O r calculations, the 
eight April 1 9S8- K ihei rawinsonde flights and the ten July 1988-Kihei launches all under 
estimate the direct r 0 and d 0 measurements by a factor of ten. (See Appendix A) This 
is expected considering the trends in the potential temperature differences. As demon- 
strated in Tigure 3, these trends generate differences approaching 2 K near the boundary 
layer top. In the stratosphere (above 15 km), potential temperature differences arc even 
larger. For this particular profile, the stratospheric potential temperature differences 
reach as much as -1 K. When these large differences arc blindly used to calculate O t . the 
result is a deceptively large over estimation of O t . The over estimated O t subsequently 
leads to an under estimation of r c and 0 0 . 

2. Case II: Using a Running Mean 

Case II investigates the effects of detrending the potential temperature profile 
prior to calculating O r The two representative high pass filters selected for this case 
include the equally weighted 1 1 1 Running Mean (RM) 
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(17) 



A) = 



lA;., + \Xj + 1 A’ + , 



and the weighted 13531 Running Mean 



X- 



1 A' ; — 2 + 3A'-_, + 5 A 7 + 3 X l+l 4- lA' +2 
13 



(18) 



Figure 4 displays the potential temperature difference profile after the 1 1 1 RM has been 
subtracted from the original potential temperature profile. There is still an increase in 
potential temperature differences around the boundary layer, however, the peak values 
of approximately 0.3 K is balanced with an adjacent value of about -0.3 K. Even the 
stratosphere appears nearly symmetrical between + - 0.7 K. 

a. Ill Running Mean 

The implementation of the equally weighted (111) RM handles the 
DigiCORA's two- to ten-second sampling rate transition by using values equally spaced 
in time when available. When data are unavailable, preference is given to the "central" 
point. Profile endpoints are ignored. For two-second data the filter takes into account 
layers between 15-30 meters thick. Ten-second data involves layers of 90-120 meters. 

b. 13531 Running Menu 

The 13531 running mean filter accommodates the two- to ten-second sam- 
pling rate change with preference given first to data equally spaced in time, and second, 
to the central point. Endpoints are noted, but generally ignored. 

c. r u Results 

Both RM methods produce r 0 values within the same order magnitude as 
the optical sensor values. With few exceptions, the relatively large and small magnitude 
variations present in the optical sensor data arc traced b\ both RM results. The 1353 1 
RM generate values approximately 16°o larger than the 1 11 RM. The error between the 
111 RM and sensor results ranges from 13% to 55%. The 13531 RM under estimates 
the sensor values by 4% to 48%. When the Q profile detrended by the 111 RM is 
truncated just below the stratosphere's large Q increase, /•„ increases by only 11%. Ta- 
ble 2 displays a sample of these results. 

d. 0 0 Res id is 

While large scale magnitude changes recorded by the sensor arc virtually 
ignored, the calculated 0 0 values are within the same order of magnitude as the sensor. 
The ranee of sensor measurements is from 3 to 7.5 u rad. Values calculated from the 
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0C05 UTC 20 April 00 




Potential Temperature Difference (Ii) 



Fig. 3. Profile of Potential Temperature Differences 
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0605 UTC 21 3 April 60 




Potential Temperature Difference (I<) 



fig. 4. I'otcn( ini Temperntuie Differences Detrended by I he III Running Mean 
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Ill RM cover a range of 0.7- 1.3 /j rad. (The 8.5 fx rad of flight 5, July 8S is ignored.) 
0 O calculated with the 13531 RM ranges from 0.8- 1.5 q rad. Truncating the detrended 
Cy profile just below the stratosphere's large Cy increase, increases 8 0 by 78%. A sample 
of these results is in "fable 2 . 



Table 2. SINGLE PROBE RESULTS - 0605 UTC 28 APRIL 1988 



C 2 r CALCULATION 
METHOD- 


r 0 (mm) 


0 0 (// rad) 


SPTD (all values) 


8 


0.2 


BG (all values) 


7 


0.2 


1 1 1 RM (all values) 


36 


1.3 


111 RM 

(troposphere only) 


40 


2.3 


13531 RM (all values) 


3S 


1.5 


SENSOR VALUES 


43 
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3. Case III: Using the Bulk Gradient Cy Equation 

The third approach incorporates an expression developed by Tatarski (1961) 
that relates Cy to bulk parameters of the atmosphere. This expression which will here- 
after be referenced as the "Bulk Gradient C 2 r " equation (BG) is. 

C 2 T =ay.Lo ; \dOlch)\ (19) 

where a is a constant: a is the ratio of eddy thermal diffusivity to eddy viscosity; L 0 is the 
outerscale length; and. (dQ'dz)- is the potential temperature gradient. Because L 0 is not 
universally well-defined and dQjd: is also vague with respect to the resolution needed to 
properly define the derivative, the values calculated from BG are suspect. The main 
strength of the expression is that it accentuates the strong dependence C\ has on the 
square of the temperature gradient. This term also provides the "atmospheric density 
inhomogcncity susceptible to being mixed by the wind turbulence" (Walters and Kunkcl. 
19S1 ). Assumptions made when applying this equation include the following: 

1. The atmosphere is mechanically turbulent. Note the lack of a shear term. 

2. The outerscale. L 0 , is half the distance between samples. 

3. Constant, a, is 2.8. 

4. Ratio, a. . is 1 . 
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a. Results 



The BG r 0 and 6 0 results are one-tenth the sensor measurements. Both op- 
tical parameters calculated in this case closely trace those values generated by the Simple 
Potential Temperature Differences (SPTD) Case. r 0 magnitudes van over a range of 6-8 
min. 6 0 values cover a range of 2-3 mm. These values are within 7% of the SPTD case. 
Considering the dominant term of the BG equation is the potential temperature gradient, 
it is not surprising to discover the final results are so closely associated to Case I. An 
explanation for why r 0 and 6 0 values are off by a factor of 10 begins with a re-evaluation 
of the assumptions. The constant, a, has a suggested range between 2.8 to 3.2 . Shifting 
this to the higher constant will not accommodate the discrepancy. Changing the 
outerscale length to the full thickness of the layer serves only to decrease r 0 and 6 0 . 
Justifying an outerscale smaller than half the thickness of the layer sampled taps into the 
fundamental problem with the Q calculations. That is, vertical resolution. The sampled 
values are too far apart to provide the temperature gradient information needed to drive 
this C 2 r equation properly. Decreasing the outerscale length without simultaneously de- 
creasing the thickness of the sampled layer is illogical. As in Case I, the BG results in- 
dicate a need for either detrending the data prior to calculating O t or using a rawinsonde 
dataset capable of producing more detailed temperature gradient information. 

4. Case IV: Using Richardson Number 

Cj- is a constant representing the total amount of energy in a turbulent layer 
(Clifford. 197S). The vertical temperature samples measured by the rawinsonde are not 
always taken from turbulent layers. To computationally isolate only those atmospheric 
layers which are turbulent, the Richardson number. Ri. was used. Richardson number 
is "the ratio of work done against gravitational stability to energy transferred from mean 
to turbulent motion" (lluschke, 197U). When Ri < 0.25, the layer is unstable and Q- is 
calculated using the SPTD. Running Mean and BC> approaches, respectively. When Ri 
> 2. the layer is stable and no C 2 r calculation is made. "Critical" conditions exist when 
0.25 < Ri < 2.0 . As the final case reviewed. O r is calculated when Ri < 2.0 . 
a. Results: Unstable Ri Case ( Ri < 0.25 ) 

Letting a turbulent layer be defined as one having a Ri < 0.25, each of the 
above three methods for calculating O t was conditionally computed. With remarkable 
consistency, all three methods over estimate ; 0 . r 0 values farthest from the sensor's value 
generally are produced by the SPTD case. Except for one extreme case where /•„ is cal- 
culated as gi eater than 360 mm from the sensor value, the SPTD r 0 values are generally 
2 to 3 times too large. The BG method leads to values that are approximately 7% closer 
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than SPTD to the sensor value. The R.M method(s) r 0 values over estimate the sensor 
by IS to 71 mm for the 13531 R.M and S to 5S mm for the 1 11 RM. 

0 0 calculations generate less consistent results. The sequence of greater to 
smaller magnitudes seen in the r 0 results does not hold for 0 0 . About 40% of all the 0 0 
data calculated for Case IV ( Ri < 0.25 > over estimate the sensor's measurement; 50% 
under estimate 6 0 ; and, 10% straddle the sensor value (see Appendix A) . 
b. Results: Conditionally Unstable Ri Case (Ri < 2.0) 

With the inclusion of more atmospheric layers in the r 0 integral, one would 
expect the resulting magnitude to be smaller. For the most part this is true. The SPTD 
and the BG cases respond in the same manner. Their r 0 values under estimate the sensor 
by 17-S0% for the SPTD case and 22-65% for the BG case. The 111 and 13531 RM 
calculations still over estimate r 0 with the 13531 R.M slightly higher than the 111 RM 
(Sec Appendix A). The magnitude of the results as compared to the sensor's value is 
generally within 78% for 13531 R.M and 59% for the 111 RM. 

0 0 results display a pattern consistent with the "All Values" results. The 
largest to smallest magnitudes calculated arc: (Sensor Value), 13531 R.M, 111 R.M, 
SPTD and BG. All values are less than 3 rad! As before the large magnitude vari- 
ations observed in the sensor data are virtually ignored in the calculated data. 

C. CORRELATION OF SENSOR AND CALCULATED PARAMETERS 

Correlating the measured and calculated optical parameters for each case discussed 
in section B generates r-squared values between 0.0009 and 0.5 . Table 3 lists r-squared 
values for each case, Note: r-squared is 1 for a perfect correlation. 

The lack of good correlation demands a review of the sensor's ability to sample 
values for Cf calculations. After detrending the potential temperature profile with a 1 1 1 
R.M. potential temperature differences for the two-second data are plotted. (Sec 
Figure 5) Initially, this profile appears to contain valid turbulence information. Re- 
peating this procedure with the temperature (K) profile (Figure 6). however, unveils a 
major problem. The systematic + -0.03 K oscillation dominating the profile implies 
that the temperature differences below 2 km are predominantly quantization noise from 
the 0.1 °C temperature resolution. Comparing Figure 5 and Figure 6shows that the 
potential temperature calculation introduces additional errors that obscure the quantiz- 
ing error. Rather than pursue a detailed correlation analysis, attention is diverted to 1) 
identifying and quantizing the major sources of error, and 2) providing recommendations 
for a temperature sensor system more suited for calculating Cf . 
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fig. 5. Initial 2 km of Detrended Potential Temperature Differences 
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Table 3. MEASURED AND CALCULATED PARAMETERS CORRELATED 



R-SQUARED 

VALUES 


r 0 (nun) 


0 o {jj. rad) 


C\ CALCU- 
LATION 
METHOD 


ALL 

Values 


Ri <0.25 


Ri < 2.0 


ALL 

Values 


Ri < 0.25 


Ri < 2.0 


SPTD 


0.009 


0.3 


0.5 


0.05 


0.3 


0.02 


BG 


0.009 


0.3 


0.4 


0.05 


0.3 


0.02 


111 RM 


0.005 


0.0009 


0.06 


0.01 


0.2 


0.1 


13531 RM 


0.002 


0.001 


0.07 


0.07 


0.2 


0.2 



D. ERROR ANALYSIS 

1. Temperature Probe Time Constant Correction 

In both the April SS and July SS soundings, the temperature was measured with 
a capacitive bead. For sensor specifications see Chapter 4. To calculate the sensor re- 
sponse time r . the following equation was used: 

1/t = A 4- Byfv , (20) 

where v is velocity, A is the constant for the calm atmosphere (r = 0) and B is a constant 
at velocity, r . 1 or this particular probe, the sensor time constant at a velocity of 6 m s 
is 2.5 s. With an ascent rate of 4 m s, t increases to approximately 3 s. In order to 
calculate the error induced by the temperature sensor time constant, the structure func- 
tion of the probe system with limiting scale lengths must be compared to the structure 
function over all scale lengths. 

Mathematically, the structure function over all scale lengths is (Tatarski, 1961): 



D(r.0,oo) = 4 [* sm : (kr/2)k~ 5/3 dk, (21) 

*'o 

where k is 2 - 7. . 7. the turbulent wavelength and r the probe separation distance. 

The structure function filtered by a one pole, low pass filter is (lluntcn. 1964): 



D{r,a,b) — 4 



( 22 ) 



r %inV/;m - s ' 3 
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dk. 



where 

a is 2 nlL 0 , the limiting lower frequency (outerscale, L 0 ), 

b is 2: r// 0 , the limiting upper frequency (inner scale, / 0 ), 

co is Ah' where V is the wind velocity (rawinsonde's ascent rate) and 

r is the time constant of the temperature probe (the upper frequency cutofi). 

The cor term in the denominator attenuates the high spatial frequencies. Since the A -5 ' 3 
term increases at low wavenumber, Eq. (22) is much more sensitive to the selection of 
an outerscale than the innerscale. Test calculations of the D(Sm.a.b) integral confirm 
that an innerscale of 0.1 m verses 0.001 nr produce virtually identical results. When 
varying the outerscale from 10 m to 1000 m (all other parameters of D(Sm.a.b) remain- 
ing unchanged), the 1000 m structure function is 200 times larger than for the 10 m. 

With an inner scale of 0.1 nr. an outerscale of 1000 nr and a 3-second temper- 
ature probe time constant, the DigiCORA's 2-second data is 15% of the actual Q 
present. The structure parameter Q, calculated from the 10-second data is 54% of the 
actual O r value. While the 10-second data appears to produce results more than three 
times closer to the desired values, this is somewhat deceptive. The larger time interval 
between samples (accent rate is constant) results in a greater probe separation distance. 
Consequently, higher spatial frequencies or rapid changes in the potential temperature 
will be missed. This aliasing of temperature measurements is one of the major difficulties 
with calculating structure function parameters from rawinsonde data. 

Another problem with calculating a structure function D from 10 second (r=40 
m) data involves a limitation of the Kolmogorov C 3 slope model. D is a function of the 
probe separation distance, r. and the outerscale. For r values comparable to the 
outerscale. the actual D(r) deviates from Kolmogorov's model. When r is larger than 
the outerscale, the measured (actual) structure function will be smaller than the D cal- 
culated from Kolmogorov's model. (See Figure 7) The significance of this problem is 
that the calculated D under estimates the actual turbulence present when r is greater 
than the outerscale. Thus, accepting the calculated Q as 54% of the actual Cf present 
is done cautiously. 
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I ig. 7. Kolmogorov Model vs. Actual Structure Function Curve. 
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a. Application of the Filter 

Correcting the basic O t equations for the temperature probe time constant 
produces an even greater over estimation of atmospheric turbulence. The magnitudes 
of the filtered r 0 and 6 0 (all values) decrease 35-50% from the unfiltered values. 
Table 4 and Table 5 display the results taken from two ravvinsonde datasets: 0605 UTC 
28 April 1988 and 1434 UTC 30 April 1988. Imposing this correction onto the calcu- 
lation of the optical parameters where Ri<0.25 and Ri<2.0 produce the same 35-50% 
decrease as the "All Values" cases. Unless otherwise noted, further reference of cor- 
rection implementation will be done on the 0605 UTC 2S April 19S8 Kihei, Maui 
rawinsondc dataset. 



Table 4. CORRECTING FOR TEMPERATURE PROBE TIME CONSTANT 



0605 UTC 28 April 198S 
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BG-A11 Values 


7 


5 
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0.1 
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36 


20 


1.3 


0.9 


13531 RM-A11 Values 


38 


22 
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SENSOR VALUES 


4S 


48 
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Table 5. CORRECTING FOR TEMPERATURE PROBE TIME CONSTANT 



1434 UTC 30 April 1988 


)'n (mm) 


O 0 [p rad) 


O t Calculation Method 
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rected 


Corrected 


Uncor- 

rected 
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SPTD-A11 Values 
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BG-A11 Values 


S 
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0.2 


0.1 


1 1 1 RM-A11 Values 


> 


16 


1.2 


0.8 


13531 RM-A11 Values 


J / 


19 


1.3 


0.9 


OPTICAL SENSOR 
VALUES 


54 


54 


7 


7 
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2. Sensor/ Atmosphere Interaction Noise 

Since the optical parameters calculated by the rawinsonde are smaller than the 
optical sensor values and compensation for the probe response-time accentuates this 
problem, some other large error source must be present. Noise introduced by the 
sensor-atmosphere interaction is such a source. Two methods attempt to quantize this 
phenomena. 

a. Standard Error Formula 

The first approach tackles the problem by focusing on the limitations of the 
rawinsonde system. In other words, knowing the resolution of each variable measured, 
how well can the given instruments measure atmospheric turbulence? This method uses 
the standard error formula which for a function E(x,y,z) is (Bevington, 1969): 
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Applying the standard error formula to the SPTD calculation of Cf produces 
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Note: T 2 and 7, are assumed to generate independent random errors. 
With Of defined as 
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the svstem induced error for Cl is 
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(26) 



Implementing the error correction reveals three layers where (Of — a c2 ) < 0. These re- 
sults imply that the measurement uncertainties are comparable to the measurements. 
1 he sensitive regions include 0 to 2 km. 4.5 to 9.5 km and 12 to 14 km. However, the 
area of major concern is near the surface where the majority of turbulence is measured. 
Eliminating the 0-2 km layer would result in a gross under estimation of turbulence 
which ultimately generates large r 0 and 0 O . Since by definition Of can not be less than 
zero, all negative error corrections are assumed to be zero. 
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The results from the Q error calculation show values less than or equal to 
0 near the surface (0 - 2 km). Between 3.5 - 9.5 km. as well as. 10 - 14 km, the Q < 0 
condition is intermittent. Unlike the Cy calculation, Q has a pressure term that effec- 
tively minimizes the influence of any values above the tropopausc. Thus, the main 
concern is with turbulence measurements near the surface. 

The average Cy structure parameter, including sensor error (0605 UTC 28 
April 198S), is (1.1 + l.S).tlO -3 A?m -2 '' 3 for the 0 - 2 km two-second data; the average 
Q ± sensor error over this same region is (S + 1 3).Yl0 _16 /?r 2 ' 3 . The 10-second data (2 km 
and above) produces an average O t ± sensor error of (7.4 + 0.8)xl0~ 2 A-m~ 2,3 ; the average 
Q, including sensor error for the 10-second data, is (2.0 + 0.4).vl0 _15 /?r 2 3 . Clearly, the 
0-2 km data are buried in uncertainties. 

b. 1 ertical Cy Quantization Noise Simulation 

The second technique used for quantizing the roundoff error in the SPTD 
(Simple Potential Temperature Differences) calculation of Cf simulates the atmosphere's 
temperature profile with a straight line slope of -0.007 K m. '1 his technique (Walters, 
1989; includes a correlation between adjacent temperature samples that is ignored in the 
standard error approach. The 0.1 K temperature resolution produces a quantization 
noise of 1.4„\T0~ 3 A-m* 2 3 for the two-second data and 6.S.rlO' 3 A’ 2 /?r 2 3 for the ten-second 
data. The two-second noise is 22° o less than the Standard Error Function's result. The 
ten-second value is the same. Running the simulation on a detrended potential temper- 
ature profile, the calculated noise for Cy is a factor of 10 smaller than the SPTD calcu- 
lation. Specifically, the two- and ten-second quantization noises are 6.2.vlO _4 A-//r 2 ' 3 and 
1 .4.V 1 0 -4 A~/?r -2 3 . respectively. 

3. Implementing Noise and Filter Corrections 
The measured Cy is 

Cj- — {Filter x ci ) + Noise, (27) 

J measured 1 actucl x J 

where the Filter is the temperature probe time constant correction, and Noise is the 
vertical Cy quantization noise from the simulation. 

Applying the noise correction to the 111 RM Cy calculation for the 0605 FTC 
28 April 198S rawinsonde flight, the original 36 mm r 0 calculated by this method in- 
creases to 46 mm (all values included). The sensor value is 48 mm. Replacing all neg- 
ative 0 T with zero reduces this measurement to 41 mm. Imposing the temperature probe 



33 



time constant correction reduces this value again to 36 mm (all values included) and ul- 
timately r 0 is 25 mm including only positive Q values. 

6 0 displays a much smaller range of resuits due to the equation's sensitivity for 
values at the tropopausc coupled with a lack of large variation created by the noise and 
filter at that level. Specifically. 0 o starts as 1.3 n rad (only the 111 RM applied). In- 
cluding all values and subtracting out the quantization noise increases 6 0 to 1.4 n rad. 
Removing all negative Cy docs not change the calculated measurement. Dividing the 
noise-free O t by the filter leads to a value of 1.0 n rad (all values). And finally, using 
only positive Q-. the d 0 result remains 1.0 [i rad. The sensor's measurement is 3 ^ rad. 

E. DETERMINING AN ADEQUATE TEMPERATURE RESOLUTION 

The (0.1 + 0.05) A' temperature probe resolution dominates the error term in Eq. 
(24). An important consideration is to determine the temperature resolution required for 
calculating Cy from rawinsonde data. First, let ct c ^ = LvIO^Awr 2 3 . A O r noise level of 
LvlO^Am' 2 ' 3 is adequate for most measurements. Assuming the temperature difference 
is 0.009S K over a 1 meter height difference (dry-adiabatic lapse rate), and that 
a. = cr = 0. one can solve for cy. Note: c\, = cy v The result is a\= 1.30.vl0 -‘AAir 3 3 . 
The roundoff error (noise) temperature variance expression (Bose, 19S5) 



where d is the number of base ten significant digits to the right of the decimal point re- 
maining after roundoff, permits the calculation of significant digits required for a given 
variance. Solving for d. the resolution necessary for a c cl . = l.vlO 4 A’m 2 3 is approxi- 
mately two decimal points remaining after roundoff, 0.U1K resolution. Repeating the 
above procedure with a A /'/Ac of -0.007 K lm (a lapse rate observed in the Kihei data- 
set), the results are also approximately two. 

The initial use of the standard error formula defined a T as 0.05 K. as implied by the 
sensor specifications. The roundoff error temperature variance expression calculates o T 
as 0.03 K. When 0.03 K is used in the Cy standard error formula, the product is ap- 
proximately one-third the magnitude of the standard error formula using 0.05 K. This 
decrease in the Cy noise term increases Q. A larger Cf leads to a larger Q and. subse- 
quently a smaller r 0 and 0 0 . Since the calculated optical parameters are already under 
estimating the sensor values, this adjustment only confirms that the current temperature 
sensor system is at best marginally acceptable for calculating Cy . 
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F. RUNNING MEAN FITTING ERROR 

A fitting error around sharp gradients, such as across inversions, arises when im- 
plementing the running mean into the C 2 r calculation. The running mean is unable to 
follow the trends. Figure S displays the large difference that occurs between the running 
mean and a given profile. Though variances are used in the Calculation, the variances 
are artificially enlarged by the poor fit across the inversion. The SPTD experiences a 
form of this fitting problem that results in aliasing. (See Section D1 of this chapter). 
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Pressure -Altitude (km) 



0G05 UTC 20 April 00 




1 1 1 Running Mean Across an Inversion 






I ig. 8. 



VI. DISCUSSION AND CONCLUSIONS 



A. CORRELATED OPTICAL DATA AND METEOROLOGICAL EVENTS 

Optical and meteorological data for the Relay Mirror Experiment (RME) were col- 
lected between 19S5-19SS in Kihei, Maui, Hawaii. Reviewing the composite Kihei op- 
tical dataset reveals the Maui site as "mediocre" to "poor". While upward and downward 
trends occur, the predominant transverse coherence length values generally stay within 
the 50 to 70 mm range. 

The lack of extremely good values (e.g., 200-400 mm) is tied to the turbulent, 3-5° 
C boundary layer inversion 2-3 km above the Kihei site. Except in the winter, north- 
northeasterly trades usually cap the boundary layer. The interface between the 4 to 7 
m s trades and the boundary layer provides the mechanical turbulence detrimental to 
optimum "seeing" conditions for light propagation. The concept of the boundary layer 
introducing optical turbulence is supported by Gossard (19S5). Using fast-response 
turbulence sensors traversing a 300m tower, he shows a turbulent layer of approximately 
0.01 F?nr 2 3 to exist just below the inversion. This implies that the greatest mixing of air 
masses with different temperatures (turbulence) is at the base of the boundary layer in- 
version. 

Kihei has a land sea breeze even during the "trade wind" season (see Chapter 2. 
Section C). This daily shift of circulation generating a warm, convective air mass over 
Kihei in the daytime followed by a cool, relatively dry environment at night renders the 
mediocre »•„ measurements. Katabatic flow off of the 3.1 km high Mt. llaleakala is an- 
other turbulence producing mechanism observed at the Kihei site. 

Profiling the atmosphere over Kihei during summer conditions, one expects a 
land sea breeze circulation, a boundary layer inversion, trade winds, a second moist layer 
aloft and a moderate jet stream. The shear surrounding any jet is a source of turbulence 
exhibited in the r fl magnitudes. Shear provides the mechanical mixing necessary for op- 
tical turbulence, however, without the compliment of a strong temperature gradient, the 
effect on the optical regime is weakened. '1 his is best exemplified by reviewing the 
summer trade season optical measurements collected when an upper-level jet is travers- 
ing the site (weak temperature gradient). The 9 0 measurements, which are sensitive to 
turbulence in the 300 to 200 mb jet stream level, are rarely below 3 n rad. Note: 3 n rad 
is considered poor optical conditions. The exception to the rule for Maui occurs when 
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a synoptic-scale event crosses overhead, tightening the temperature gradient associated 
with the jet (see Appendix D). A stronger example is examining the influence of a polar 
jet. The tighter temperature gradient and subsequently stronger winds associated with 
the polar jet are able to drive 0 0 measurements down to 1 to 2 n rads, a very turbulent 
environment optically. An example of this was in the 10 December 1985 Haleakala data, 
where 0 0 was measured as 1.6 n rad. 

One of the most optically active Hawaiian cases occurs in April 1988, a winter case. 
The Kihei site experiences the passage of a low or "eddy". A detailed discussion of the 
optical and meteorological conditions throughout that experiment is found in Appendix 
D. 

B. OPTICAL VARIABLES CALCULATED FROM RAWINSONDE DATA 

During four of six RML experiments at Kihei, Maui on-site rawinsonde launches 
were released. Two sessions employed the use of a Vaisala DigiCORA MW 11 
Rawinsonde System. The eight April 19S8 launches and ten flights recorded in July 1 9S8 
began with a two-second sampling rate up to 2 km, followed by a ten-second sampling 
rate for the remainder of the flight. All 18 vertical profile sets were utilized for testing 
four distinct Cl calculation methods. 'I he final rating of each method was done by using 
the O r profile to calculate the two optical parameters. r 0 and 0 0 . These calculated values 
were compared to the optical measurements acquired simultaneously with the 
rawinsonde launch. 

The first two methods, the Simple Potential Temperature Differences (SPTD) and 
the Bulk Gradient (BG). produced similar results. Both over estimated the turbulence 
leading to an under estimation of /•„ and 0 U by a factor of 10. 

Using the entire profile, the 1 1 1 and 13531 running mean techniques produced r 0 and 
0 0 values of the same order magnitude as the sensor values. The r 0 data even showed a 
general ability to match the relatively large and small magnitude variations observed in 
the sensor data. 0 0 calculated from the running mean methods were generally insensitive 
to the magnitude changes measured by the isoplanometer. However, the calculated 
values were within the same order of magnitude as the sensor. Overall, the 13531 run- 
ning mean provided turbulence parameters at least 16% closer to the sensor values than 
the 1 1 1 running mean. 

Using the Richardson number to isolate the turbulent layer from the noil-turbulent 
brought mixed results. Conditioning the Cf profile by calculating turbulence for layers 
with Ri < 0.25 generated /■„ two to three times larger than the sensor. The BG Case 



provided values approximately 7% closer to the sensor. The running mean approaches 
were within 8-71 mm of the sensor. 0 o calculations for SPTD. BG and running mean 
generated results such that about 40% over estimated the sensor's 0 0 ; 50% under esti- 
mated; and, 10% straddled the sensor value. 

The conditionally unstable results (Ri<2) for the SPTD Case under estimated the 
r 0 sensor by 17-80%. 

Poor correlation between the calculated and measured optical parameters revealed 
the existence of serious problems. Limitations in the current sensor system prevented 
successful Q- calculations. The limitations were broken down into three areas: 

1. Excessive Temperature Probe Response Time, 

2. Limited Vertical Resolution (aliasing) and 

3. Limited Temperature Resolution (Roundoff Error). 

The limited temporal resolution of the DigiCORA MW 11 rawinsonde profiles was 
confirmed by calculating the temperature probe time constant correction. Lor the 
2-second data, only 15% of the actual O r was measured. While the 10-second filter 
displayed 54% of the actual Cf. this was deceptive. The distance between 10-second 
samples produced aliasing effects when the mean potential temperature gradient was 
large. This inability to resolve high frequencies artificially enlarged the already over es- 
timated turbulence present. In addition, this large 40-45 m separation approached the 
Kolmogorov outerscale. 

The 2-second data was not as susceptible to aliasing. Instead, the 0.1 K temperature 
resolution (roundoff error) introduced quantization noise. Using two distinct proce- 
dures. the standard error formula and a vertical Cf quantization noise simulation model, 
noise corrections of 6.2.\T0~ 4 A’h;r 2 3 (2-second data) and 1.4.vlO~ 4 A’ 2 //r 2 3 (10-second data) 
were calculated. Application of this correction confirmed that the O t values in the initial 
2 km of data were mostly quantization noise. Because most optical turbulence is ex- 
pected near the surface, this realization was not encouraging. 

Another major problem with calculating C? T from rawinsonde data was insufficient 
vertical resolution that produced a fitting problem with the running mean case. Just as 
the inability to resolve higher frequencies (aliasing) produced an over estimation of O t 
for the SPTD Case, this effect became even more apparent when the running mean 
methods encountered a 3-5 0 C inversion with a 10-second sampling rate. This enhanced 
the existing over estimation of turbulence. In contrast, turbulence calculated across the 
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sharp temperature shift (inversion) captured within the 2-second data (1537 UTC 7 July 
19SS) did not appear to suffer from this difficulty. 
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VII. RECOMMENDATIONS FOR FUTURE RESEARCH 



A. LINKING OPTICAL DATA AND METEOROLOGICAL EVENTS 

Quantizing the optical turbulence based on weather conditions is a non-trivial 
project. As simultaneous optical and meteorological observations accumulate, the pos- 
sibility of plotting "weather" charts of iso — r 0 and iso — 6 0 lines nears reality. Such a 
task, however, requires at least two phases of investigation. The first deals with 
synoptic-scale events. Major atmospheric phenomena affecting optical propagation in- 
clude: the passage of cold and warm fronts, the polar and subtropical jets, high and low 
pressure systems, and so on. A dedicated climatological study (spanning many years) 
of optical parameters during each of these events, with measurements taken at one lo- 
cation. would provide a statistical sample from which ranges of values could be deter- 
mined. Being realistic, a more cost effective and time cfTicicnt approach for acquiring a 
(marginally) valid statistical sample is to design a set of experiments with multiple 
(minimum 3. preferable 5) r 0 and 0 0 sensors. T he geometric shape formed by 4 of the 
sensor site locations would be a square. The sides of this square would be along the 
prevailing winds and measure a kilometer or more in length. The fifth sensor would be 
in the center of this box. Time based cross-sections would be derived from three of the 
sensors. Because the central site is common to both possible cross-sections, this sensor 
would be the most critical. 

The meteorological sensors needed for simultaneous measurements would depend 
on the particular event being studied. T or example, the passage of a cold front would 
require surface meteorological stations, rawinsondes and satellite resources. The effect 
of a polar jet would require rawinsondes. aircraft and satellite data. T he main objective 
is to pinpoint optical measurements (processed in real time) with the movement of the 
weather phenomena. Armed with at least three "clean" case studies of a specific weather 
event, r 0 and 0 0 ranges for various sectors of the specified event would be possible. 

The second area of investigation involves the meso- and microscale weather activ- 
ities. Though weather phenomena on these scales arc often site specific, one major 
source of optical distortions common to all sites is the boundary layer. The influence 
of the boundary layer has been estimated as 30-50%. A simultaneous data collection 
at the surface and above the boundary layer inversion would better quantize the impor- 
tance of this atmospheric structure. Because Mt. Haleakala is over 3 km high, such a 
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project could theoretically be done in Maui. However. Haleakala's Crater experiences 
a local jet. Therefore, the measurements above the boundary layer inversion would have 
to be sampled sufficiently far from this accelerated How. Another possible location is 
along the central California coast where r 0 values have often been recorded to be more 
than 200 mm. A dedicated study contrasting the good "seeing" conditions above and 
below the boundary layer interface may provide a better resolution for this question. 

Another meteorological phenomena at the mesoscale is gravity waves. With the 
current method of data acquisition for the isoplanometer, tracing the influence of gravity 
waves in real-time is possible. This ability coupled with the mobility of the 
isoplanometer could be of value to gravity wave studies. 

B. CALCULATING Cy FROM RAWINSONDE DATA 

Another approach to resolving the correlation of optical turbulence and 
meteorological conditions is to calculate key optical parameters (such as r 0 and 0 Q ) di- 
rectly from rawinsonde data. Currently, one of the main obstacles is the need for a low 
cost temperature probe and operating system capable of acquiring large high resolution 
data bases. An adequate single probe rawinsonde temperature sensor would require a 
0.0 1-0. 02 K resolution. 100 msec response time and a 1 m vertical resolution. These 
specifications are not unreasonable with current technology. Some method of detrend- 
ing the potential temperature profile with a running mean still will be necessary until 
0. 1-0.2 nr vertical resolution is possible. 

An interesting experiment would be to launch a rawinsonde with both the above 
recommended single temperature probe and the more expensive two temperature probe 
sensor. The latter has been discussed in Olmstead (19SS). A comparison of the two O t 
profiles generated from these methods would help to calibrate the less expensive single 
probe method. 

As personal computers become faster and high data-acquisition rate rawinsonde 
packages become routine, the theoretical concept of evaluating and predicting the op- 
tical turbulence o\er any site in real time will evolve and become a reality. 



APPENDIX A. CALCULATED/MEASURED OPTICAL PARAMETERS 



The optica] variables r 0 (transverse coherence length) and 0 O (isoplanatic angle) de- 
pend on an integral of Q along the vertical path. The refractive structure parameter Q 
is proportional to the temperature structure parameter O r which can be computed from 
the rawinsonde data. Four methods for computing the temperature structure parameter 
O t from the rawinsonde data include: 

1. the Simple Potential Temperature Differences (SPTD) Equation, 

2. detrending the potential temperature differences with a 111 Running Mean (111 
RM) or a 13531 Running Mean (13531 RM), 

3. the Bulk Gradient Equation (BG) and 

4. conditional data sampling based on the Richardson Number (Ri<0.25 and 
Ri < 2.0). 

Each O t method was tested by calculating /„ and 0 0 from the O t profile. These results 
were ultimately compared to the optical measurements simultaneously acquired over the 
same site. 

Table 6 and Table 7 correlates the (light numbers referenced in the Appendix A 
plots with their rawinsonde launch time (ETC). 



Table 6. APRIL 19, SS KIHEI. MAUI RAWINSONDE LAUNCHES 



FLIGHT 

NO. 


TIME (ETC) 


DATE (DAY MON YR) 


i ; 


OMO 


27 April 198S 


2 


1548 


27 April 1988 


■s 


0605 


2S April 1988 


4 


1549 


2S April 19S8 


5 


0005 


29 April 19SS 


6 


155S 


29 April 19SS 


7 


0544 


30 April 1988 


s 


1434 


30 April 19SS 
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Table 7. JULY 1988 KIHEL MAUI RAWINSONDE LAUNCHES 



FLIGHT 

NO. 


TIME (UTC) 


DATE (DAY MON YR) 


1 


0703 


7 July 198S 


2 


1537 


7 July 19SS 


-> 

j 


0536 


S July 19SS 


4 


1612 


S July 19SS 


5 


0605 


9 July 19SS 


6 


1527 


9 July 19SS 


7 


0556 


10 July 19S8 


S 


1532 


10 July 19SS 


9 


0550 


11 July 19SS 


10 


0926 


11 July 19SS 



4 - 



A. R<, RESULTS 

Figure 9 through Figure 1 1 display the /•„ values calculated from the eight April 
1988 and ten July 198S Maui rawinsonde datasets. An estimate of the r 0 simultaneously 
measured over the same Kihei. Maui site is also included. 

For the measured r 0 time series corresponding to these profiles, sec Nelson (198S). 
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Fig. 10. Mcnsuied and Calculated r 0 (Hi < 0.25) 
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fig. 11. Measured and Calculated r„ (Ri < 2.0) 
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B. 0 0 RESULTS 

Tigure 12 through Figure 14 display the 6 0 values calculated from the eight April 
1988 and ten July 1988 Maui rawinsondc datasets. An estimate of the 0 0 simultaneously 
measured over the same Kihei, Maui site is also included. 

For the measured 0 D time series corresponding to these profiles, see Nelson (1988). 
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f ig. 12. Measured and Calculated 0 0 (All Values Included) 
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fig. 13. Measured and Calculated 0 0 (Ri < 0.25) 
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Fig. 14. Measured and Calculated 0 0 (Ri < 2.0) 




APPENDIX B. MAUI, HAWAII - AUGUST 1987 



A. METEOROLOGICAL DATA 

Appendix B includes the August 1987 vertical profiles of atmospheric data acquired 
during the Relay Mirror Experiment at Kihei, Maui. The data acquisition system used 
was a VIZ WO-SOOO RP+ rawinsonde system. Vertical resolution is approximately 35 
m. Launch site was at 20.45° latitude, 156.27° longitude. Table 8 couples the 
rawinsonde flight numbers and launch times for all the August 19S7 Kihei soundings. 



Table 8. AUGUST 1987 KIHEI. MAUI RAWINSONDE LAUNCHES 



FLIGHT 

NO. 


TIME (ETC) 


DATE (DAY MON YR) 


1 


0658 


18 August 1987 




1440 


IS August 1987 




0722 


19 August 1987 


4 


1103 


19 August 1987 


5 


1512 


19 August 1987 


6 


0607 


20 August 19S7 


7 


loll 


20 August 1 9S” 


8 


0650 


21 August 19S7 


9 


1106 


21 August 1987 


10 


1559 


21 August 1987 


1 1 


0541 


22 August 1987 


12 


105? 


22 August 1987 
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fig. 15. Kiliei. Maui Rnwinsoude Data: 0658 UIC 18 August 1987 
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August 07 - Flight 2 
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Fig. 16. Kiliei. Maui Rnuinsoncle Data: 1440 UTC 18 August 1987 
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Tig. 17. Kiliei, Maui Rauinsmule Data: 0722 U IC 19 August 1987 
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Fig. 18. Kiliei. Maui Rnuinsondc Data: 1103 UTC 19 August 1987 



57 



August 07 - Flight 5 
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fig. 19. Kilici. Maui Rauinsomle Data: 1512 UTC 19 August 1987 
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Kiliei. Maui Raw insonde Data: 0607 UTC 20 August 1987 



59 



August 07 - Flight 7 



Mixing Ratio ( g/kg ) 




Wind Direction (deg) 




Fig. 21. Kiliei. Maui Rauinsonde Data: Mill UIC 20 August 1987 
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Fig. 22 . Kihci, Maui Rnuinsonde Data: 0650 IJ'TC 2 \ August I9<S7 
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Tig. 23. kihei. Maui Rawinsonde Data: 1 1 DO UIC 21 August I9S7 
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fig. 24. Kiliei. Maui Rauinsonde Data: 1559 UIC 21 August 1987 
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August 07 - Flight 11 
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Fig. 25. Kihei, Maui RaMnsonde Data: 054 1 UTC 22 August 1987 
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August 07 - Flight 12 
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fig. 26. Kiliei, Mntii Raninsonde Data: 1053 UTC 22 August 1987 
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APPENDIX C. MAUI, HAWAII - DECEMBER 1987 



A. METEOROLOGICAL DATA 

Appendix C includes vertical profiles from the December 1987 Kihei, Maui 
rawinsonde launches. The VIZ WO-SOOO RP+ rawinsonde system used to acquire this 
data sampled values approximately every 35 m. Launch site for these profiles was at 
20.54° latitude, 156.27° longitude. Table 9 tabulates all the flight numbers and launch 
times for the December 1987 Kihei soundings. 



Table 9. DECEMBER 1987 KIHEI, MAUI RAWINSONDE LAUNCHES 



FLIGHT 

NO. 


TIME (UTC) 


DATE (DAY MON YR) 


1 


0950 


5 December 1987 




1552 


5 December 19S7 




0639 


6 December 198" 


1 

-i 


1131 


6 December 1987 


N 


0517 


7 December 198“ 


6 


No Data 


-- 


7 


1624 


7 December 1987 


8 


0537 


S December 1987 


9 


1031 


8 December 198“ 


10 


1557 


8 December 198“ 


1 1 


0345 


9 December 198“ 


12 


095 S 


9 December 1 Q 8“ 
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Fig. 27. Kiliei, Maui Rnuinsowlc Data: 0950 IJTC 5 December 1987 
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Kilici. Maui Raninsonde Data: 0639 UIC 6 December 1987 
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fig. 30. Kilici. Maui Rnuiusonde Data: 1131 UTC f> December 1987 
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Fig. 31. Kiliei, Maui Rauinsoiule Dala: 0517 UIC 7 December 1987 
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Fig. 32. Kihei, JNlnui Rauinsonde Data: 1624 UTC 7 December 1987 
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Fig. 33. Kilici, Maui RaMnsonde Data: 0537 UTC 8 December 1987 
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December 07 — Flight 9 
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I i”. 3-1. Kihci, Maui Rauinsoude Data: 1031 IJI C 8 December I ‘>87 
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Fig. 35. Kihci, Maui RaMnsonde Data: 1557 UTC 8 December 1987 
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Tig. 36. kiliei, IMaui Raninsomle Data: 0345 UI C 9 December 1987 
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Fig. 37. Kilici, Maui Raninsonde Data: 0958 UIC 9 December 1987 
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APPENDIX D. MAUI, HAWAII - APRIL 1988 



A. METEOROLOGICAL DATA 

Appendix D includes the April 19SS vertical profiles of atmospheric data acquired 
for the Relay Mirror Experiment at Kihci, Maui. The data acquisition system used was 
a Vaisala DigiCORA MW 11 rawinsonde system. Launch site was at 20.75° latitude, 
156.43° longitude. The initial 2 km were sampled every two seconds (8-10 m). Data 
acquisition above 2 km was once every ten seconds (40-45 m). Flight number and 
launch times for all the April 19SS soundings are provided in Table 10 . 



Table 10. APRIL 1988 KIHEI. MAUI RAWINSONDE LAUNCHES 



FLIGHT 

NO. 


TIME (UTC) 


DATE (DAY 


MON YR) 


1 


0610 


27 April 


19SS 


*) 


1548 


27 April 


1988 




0605 


2S April 


19S8 


— r 


1549 


28 April 


198S 


s 


06<)5 


29 April 


1988 


6 


1558 


29 April 


198S 




0544 


30 April 


1988 


8 


1434 


30 April 


19SS 



7S 




Temperature & Dewpoint (C) 



Fig. 38. Kilici. Maui Rauinsonde Data: 0610 IJIC 27 April 1988 
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Fig. 39. Kiliei, Maui Rawiusoiide Data: 0610 UIC 27 April 198” 
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Fig. 40. Kihci, Maui Ran insonde Data: 0610 U IC 27 April 1988 
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f ig. 41. Kilici, Maui Rauinsoiule Da(n: 1548 IJTC 27 April 1988 
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42. Kiliei. Maui Rauinsonde Data: 15-18 UTC 27 April 1988 
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fig. 43. Kihci. Maui Rnuiiisomle Data: 1548 UTC 27 April 1988 
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Fig. 45. Kilici, Maui Rawinsomle Data: 0605 IJIC 28 A pi il 1988 
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Fig. 46. Kilici, Maui Rauinsonde Data: 0605 UTC 28 April 1988 
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fig. 47. Kilici, Maui Rauinsondc Data: 1549 U1C 28 Apiil 1988 
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Fig. 48. Kiliei, Maui Raninsonde Data: 1549 UTC 28 April 1988 
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Fig. 50. Kihei, Maui Rawinsonde Data: 0605 IJTC 29 April 1988 
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Tig. 51. Kiliei. Maui Rauinsonde Data: 0605 UTC 29 April 1988 
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Kiliei, Maui Rauinsonde Data: 0605 UTC 29 April 19.SS 
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g. 53. Kihei, Maui Rawinsonde Data: 1558 U IC 29 Apiil 1988 
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Fig. 54. Kihei, Maui Rawinsonde Data: 1558 UTC 29 April 1988 
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Fig. 55. kiliei, Maui Rauinsomle Data: 1558 U I C 29 April 1988 
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Fig. 56. Kihei, Maui Rawinsomle Data: 05-44 UIC 30 April 1988 
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F ig. 57. Kiliei. Maui Rauinsoiule Data: 05-4-1 IJTC 30 April 1988 
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Kilici, Maui Rnninsonde Data: 0544 UTC 30 April 1988 
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fig. 59. Kilici, Maui Rauinsonde Data: 1*134 UIC 30 April 1988 
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Fig. 60. Kiliei, Maui Ran insonde Data: 1434 UK' 30 Apiil 19‘>‘? 
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F ig. 61. Kiliei, Maui RaMnsomle Data: 143-1 UTC 30 April 19S8 
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B. APRIL 1988, EXPERIMENT REPORT 

The following is a general meteorological summary corresponding to the optical 
data measured during the 27 April 19SS (UTC) through 01 May 1988 (UTC) Kihci, 
Maui, Hawaii Experiment. 

1. 0600-1730 UTC 27 April 1988 

Extensive cloud cover prohibited optical measurements until approximately 
1100 L”I C. The 0610 UTC Kihei rawinsonde launch includes the following surface pa- 
rameters: 

P (Pressure) = 1005.9 mb 
T (Temperature) = 22.6° C 
RH (Relative Humidity) = 71% 

Wind Speed. Wind Direction = calm ( < 2m 's), N A 
Stability = unstable between surface and 0.2 km 

The stratus deck restricting optical measurements occurs between 2. 2-2. 7 km. and ex- 
tends horizontally to the NE 3-5° latitude according to the GOES- WEST Infrared (IR) 
Satellite Image taken within 30 minutes of balloon launch time. A temperature inversion 
of 5.1° C between 2.9 and 3.2 km caps the cloud layer, and is coupled with light winds 
( < 5 m s) from the NE. 

Winds gradually increase with height until 9.9 km (289.7 mb) where a maximum 
wind speed of 22.0 m s is reported. Wind direction remains from the eastern quadrants 
(0-180° from the North) up to 13.2 km (173.1 mb). Within the next 97m (174.4-171.1 
mb) moderate winds become calm and veer 201 degrees, 'flic subsequent increase in 
wind speed peaks at 13.7 m s (143.2 mb) before it gradually decreases to winds of less 
than 10 m s. Wind direction continues to be from the NW until IS. 7 km (well after the 
initial signature of the tropopause) where it crosses into the eastern quadrants for the 
remainder of the vertical sampling. The rawinsonde flight is terminated after 90 minutes 
10 seconds at 25.019 km (24.9 mb). 

Comparing the 0546 UTC and the 1246 UTC IR Satellite Images, the dense 
stratus deck moves from NE of the site to SW. The 1246 UTC and 1446 UTC IR Sat- 
ellite Images display very little change in the intermittent stratus deck, though r 0 meas- 
urements indicate a 1 cm increase over the same time interval. 

Shortly before sunrise (154S UTC 27 April 19S8) a second rawinsonde launch 

records the following surface parameters: 

P (Pressure) = 1005.3 mb 
T (Temperature) = 1S.S° C 
RII (Relative Humidity) = 87% 
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Wind Speed, Wind Direction = calm ( < 2m s), X'A 
Stability = stable conditions 

The first of two inversions is found between the surface (0.083 km) and 0.257 km. A 
second slight temperature increase of 0.4° C (0.365-0.430 km) is capped by a small hu- 
midity increase of 13% between 0.459-1.012 km. The small temperature inversion is 
supported by 0.9 m s increase in wind speed. Wind direction remains constant at 28° . 
Another increase in moisture content is displayed between 1.270-1.737 km. The RH 
maximum of 83% is reported at 1.737 km. Winds sampled thus far are light and from 
the NE, WE. For a thin layer, 1.7S2-2.275 km, the light winds back to northwesterlies. 
A sharp decrease of 52% RI1 (2.075-2.713 km) overlaps the variable winds. Embedded 
within this layer is a 1.5° C inversion between 2.511-2.824 km. Except for a slight 
moisture increase between 7.754-S.824 km. the remainder of the sounding displays an 
extremely dry atmosphere. At 9.520 km, winds reach a peak velocity of 17.7 m s. At a 
slightly higher level, 10.558-12.396 km. wind speed decreases 11.9 m s and veers 163 de- 
grees. The veering continues another 54 degrees over the next 0.943 km. The temper- 
ature gradient associated with the 217 degree wind shift is -20° C 2.781 km (2.61° C'km 
less than the dry-adiabatic lapse rate). The tropopause begins around 17.3 km. 

2. 0600-1730 UTC 28 April 1988 

Kihei rawinsonde surface measurements recorded at 0605 UTC include: 

P (Pressure) = 1005.3 mb 
T (Temperature) = 22.6' C 
RH (Relative Humidity) = 72% 

Wind Speed. Wind Direction = calm ( < 2m s). N A 
Stability = unstable between surface and 0.2 km 

RH reaches a maximum of 91 To just prior to the first of two inversions. The initial in- 
version is between 2.254-2.430 km (778.4-762.0 mb). This 2.3° C increase is coupled with 
an S degree veering of winds, and no change in wind speed. The second inversion is 
fragmented into two parts. The total layer displays a l.l&degee C increase over 0.4 km 
(2.793-3.193 km). Winds associated with this stable layer increase by 1 m s and veer 15 
degrees. As before, the inversion is preceded by a slight increase (+ 1S%) in RH. 

Between 5.627 km and 6.004 km light winds veer 132 degrees NEE to SW. The 
temperature gradient associated with this large wind shift is -3C .377 km (1.S42 C, km 
less than the dry-adiabatic lapse rate). 
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Another large veering in wind direction is between 8.951-10.830 km. This 203 
degree change is coupled with a 9.9 m s increase in wind speed and a temperature gra- 
dient of -13.2 1.S79 km (2.775 C/km less than the dry-adiabatic alpse rate). 

The jet stream occurs between 1 1.561-14.431 km and flows from the NNW.NW. 
Peak winds of 42.1 m's at 315° are evidient at 12.493 km. The temperature gradient 
across the jet is -24.7° 02.870 km (-8.606° C/km). The Lihue rawinsonde data recorded 
about 5.25 hours after the local Kihei launch displays the jet stream between 
12. 354-14. 5S0 km (wind direction across the jet backs from 292° to 279° ). Temperature 
gradient associated with Lihue's jet stream measurements is approximately -19.0° 
C 2.652 km (7.164° C km). Peak wind is reported as 42. S m s, 286° at 13.4S7 km; 42.8 
m s, 2S3° at 13.940 km (Note: r 0 is 5 cm during both sampling intervals 0605 and 1115 
UTC.) 

At 1549 UTC 2S April 19SS. surface measurements are: 

P (Pressure) = 1005.3 mb 
T (Temperature) = 17.9° C 
RI1 (Relative Humidity) = 90% 

Wind Speed. Wind Direction = calm ( < 2m s). N A 
Stability = stable 

A 2.3° C inversion, as well as a 14% drop in RII. marks the first 84 m of the vertical 
profiles. Winds remain light and variable for the initial 2 km. A second minor inversion 
is between 2. 368 and 2.S07 km displaying a temperature gradient of 0.5° C 0.438 km 
(1.134 ° C/km). The easterly winds associated with this event increase 2.6 m s and veer 
6 degrees. 

A general wind direction shift from eastern quadrants to western quadrants 
takes place around 500 mb (6.137 km). Note: 500 mb is often labelled a "level of non- 
divergence". The Level of Non-Divergence is defined as a "mid-tropospheric surface 
separating major regions of horizontal convergence and divergent associated with the 
typical vertical structure of the migratory cyclonic scale weather systems" (Iluschke, 
1970). In the upper half of the profile (above 500mb), RII increases at a rate of 36% 
in 2.664 km. This secondary peak of 49% Rll is measured at 12.005 km. Nestled above 
the humidity influx is the jet stream. 12.765-14.448 km. Peak winds of 47.0 m s are re- 
corded at 13.319-13.361 km. The temperature gradient across tire jet stream is -12.3° 
C 1.6S3 km (-7.308° C km). 
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3. 0600-1730 UTC 29 April 1988 

Optical data gathering during this session was severlv restricted by cloud cover 

from a synoptical scale weather system. Unlike previous nights, the 0605 UTC sounding 

reveals the near surface winds to be from the S\V. (This implies the presence of an 

anabatic flow between Kihei and Haleakala.) The rest of the surface parameters are: 

P (Pressure) = 1006.5 mb 
T (Temperature) = 24.0° C 
Ril (Relative Humidity) = 75% 

Wind Speed. Wind Direction = calm (< 2m's). X A 
Stability = unstable between surface and 0.2 km 

The vertical humidity profile displys two peaks. The first, a gradual 13% km 
increase climaxing at 84% between 1.261-1.330 km. A minor 0.1° C 0.00S m inversion 
caps this first humidity maximum. The second peak follows a sharp 73%. km gradient 
that sends the humidity soaring to 75% at 8.976 km. Most significant to this sampling, 
however, is the wind profile. As stated before, near surface winds (0.219-1.682 km) are 
light and SW. A backing into the eastern quadrants follows and persists between 
1.759-5.104 km. A 0.5° C inversion interrupts the atmospheric structure between 
2.017-2.211 km. Above 5.1 km. winds veer into the western quadrants where they re- 
main for the rest of the vertical sounding. The jet stream begins at 13.046 km and ex- 
tends to 14.902 km. The maximum velocity sampled is 39.7 m s (13.721-13.773 km). 
The temperature gradient accompanying the jet is - 1 1.4° C 1.856 km (-6.142° C km). 

Using the 0546 UTC 1R Satellite Image, the site appears in the area of a comma 
cloud where one would expect overrunning and warm advcction; the warm front would 
be to the south. At 1146 UTC. the site appears to straddle the cold front. T he Lihue 
sounding launched at 1115 UTC, records a surface pressure of approximately 3 mb 
greater than Kihei. after interpolating to S3 m (height of the Kihei site). Since Lihue is 
approximately 2° north and 3° west of Kihei, this post-cold front characteristic is ex- 
pected. Associated with the approximate time of cold front crossing (and all the ex- 
pected tightening of gradients), the r 0 sensor displays a 20+ mm drop (!) within the same 
time interval. At 1300 UTC, r 0 measurements indicate a 15 nun plunge in magnitude. 
The corresponding satellite image implies the passing of additional cloud lines associated 
with smaller fractions of the larger cold front during that time period. The 1446 UTC 
image displays the passing of the most extensive cloud cover associated with the cold 
front in the comma cloud. By 1546 UTC ( 12 minutes prior to Kihei sonde launch 6). the 
site is boardering the area that is normally (in more northerly latitudes) associated with 



106 



the dry tongue, "a protrusion of relatively dry air into a region of higher moisture con- 
tent" (Huschke, 1970). In this sector, one expects cold advection and sinking air. 

The 1558 UTC Kihei rawinsonde launch records the following surface parame- 
ters: 

P (Pressure) = 1006.6 mb 
T (Temperature) = 18.0° C 
RH (Relative Humidity) = 86% 

Wind Speed, Wind Direction = calm ( < 2m s), N/A 
Stability = stable 

A small inversion of 2.8° C over 57 m initials the profile. Winds remain light and vari- 
able throughout the first 5 km. Humidity measurements indicate a small peak of 72% 
at 0.561 km followed by two minor inversions. A second peak of 71% at 1.715 km 
precedes a small 0.5° C increase between 1.768-1.815 km. A third humidity maximum 
of 62% is recorded at 3.159 km with yet another maximum of 69% at 7.658 km. The 
highest humidity measurement (apart from the surface) occurs at 9.307 km (73%). with 
a second strong influx of humidity at 10.151 km (72%). This moist layer immediately 
precedes a westerly jet stream (10.964-13.037 km). Unlike the former soundings, this 
jet's maximum wind speed is 38.6 m s. The temperature gradient across the jet stream 
is -13.6° C 2.13S km (-6.361° C'km). Winds remain out of the western quadrants well 
past the initial signature of the tropopause (17 km). Not until 20 km, do the winds 
suddenly become light and variable. 

4. 0500-1730 UTC 30 April 1988 

The 0544 UTC Kihei launch surface parameters include: 

P (Pressure) = 1007.9 mb 
T (Temperature) = 23.3° C 
RH (Relative Humidity) = 79% 

Wind Speed. Wind Direction = light, W 
Stability = stable 

Surface stability is defined by a 0.5° C inversion in the first 6 m of the profile. An in- 
termittent layer of stratus begins at 2.082 km (96% RH). A second stable layer marked 
by a 1.3° C increase in temperature over 0.379 km (3.430° C/km) caps the moist layer. 
By 2.6S7 km the winds pick up and veer from the NW to SE. The second of three peaks 
in the humidity is recorded at 7.113 km (84% RH). The third humidity maximum at 
8.9S2 km (75%) immediately follows a layer of strong winds (33.8 m's) from the SW. 
A weak 0.6° C 0. 101 km (5.941° C km) inversion ties this layer's structure together. A 
slow westerly jet stream (max wind speed of 36.8k m s at 12.972-13.021 km) occurs be- 
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tween 12.716 km and 13.228 km. This is supported by a -2.2° C/0.512 km (-4.297° 
C km) temperature gradient. Despite a decrease in wind speed with height, winds remain 
W,SW up through 18 km. Around the height of the tropopause (20.5 km) winds back 
into the eastern quadrants where they remain until the sampling is terminated (23.975 
km). 

Consulting the 1115 UTC Lihue sounding, as well as the 1146 L'TC satellite 
image, the overall structure of the comma cloud is now fragmented. Low level clouds 
are intermittently dispersed around the Hawaiian Islands. A moist layer between 1-2.1 
km is implied from Lihue's sounding. Surface winds are moderate to light and from the 
east until about 6 km where they veer and increase in magnitude. The maximum velocity 
recorded, however, is 34.0 m s from the west (12.536 km). This is hardly the strength 
of the system just monitored less than 4S hours ago. The /•„ measurements reflect this 
lack of gross turbulence with an increase from 45 mm to nearly 70 mm sampled over the 
same 5-6 hour time period. 

The last Kihei rawinsonde launch of the experiment (1434 UTC) records the 

following surface measurements: 

P (Pressure) = 1007.4 mb 
T (Temperature) = 1S.0° C 
Rll (Relative Humidity) = 100% 

Wind Speed. Wind Direction = light. XL 
Stability = stable 

The vertical profile begins with a 3.4° C inversion in the first 36 m. Surface humidity 
of 100% (a measurement confirmed by condensation on the telescope lenses) decreases 
20% in the initial 0.312km. The light XL surface winds back to light W.SW within 
0.265km of the surface. Humidity increases to a peak of 83% at 0.638 km. dips about 
5%. then increases to S7% by 1.241 km. A 0.5° C increase occurs between 1.324-1.440 
km (4.310° C km), followed by another inversion of 0.6° C between 2.315-2.464 km. 
Winds display an increasing trend until approximately 4.6 km (SW). At this same level, 
humidity plunges to values of less than 10%. The westerly jet stream is nestled between 
12.060-13.196 km carrying a peak velocity of 37.9 m s betwen 12.444-12.521 km. The 
tropopause occurs around 19.1 km. 

5. 0600-1730 UTC 01 May 1988 

Early evening ; 0 data indicates intermittent cloud cover. The 0546 UTC satellite 
image displays cloud cover as a local event. The closest cyclone is XW of the site by 
approximately 10° latitude (1 100 km). The storm system that traveled through the area 
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earlier in the week has dissolved into the high pressure ridge Nil of Hawaii. The two 
available soundings (Hilo and Lihue) launched at 1100 UTC and 1115 L'TC. respec- 
tively, correspond to relatively high r 0 measurements (SO mm). Both rawinsonde 
soundings indicate an absence of surface inversion. Dewpoint depression minimums (no 
RH is recorded in these two rawinsonde datasets) occur between 0.4 and 1.6 km implying 
a moist layer within this 1.2 km region. Super-adiabatic laypsc rates are noted between 
surface- 1000 mb (12.327° C km) at Lihue. and 709-704 mb (46.867° C km) at Hilo. 
Because there is no consistent level for these extremely large lapse rates, they are most 
likely locally induced. Winds below 4 km are generally from the eastern quadrants 
(Lihue’s easterly winds extend to 5 km). Above this level, wind speed strengthens and 
wind direction is from the west. Jet stream velocities are recorded as W, NW between 
11.792-13.922 km over Hilo, and W.SW between 10.645-14.672 km over Lihue. Maxi- 
mum wind speeds of 42.8 m s at 1 1.785 km (Hilo) and 39.2 m s at 12.016 km (Lihue) 
are recorded. Both sites measure a lapse rate of -6.7° C km across the jet. The two 
soundings also concur that the tropopasue caps the jet's subsiding winds at approxi- 
mately 17.3 km. 

/•„ measurements record a gradual 30 mm decrease between 1200-1300 L'TC. 
The only available resource over that time interval is satellite images. Apart from an 
apparent increase in local cloud cover, there is no glaring evidence of turbulent sources. 
A review of 19SS May 2 reveals the development of another comma cloud structure west 
of Hawaii. Perhaps the sudden drop in / 0 (turbulence increase) is a results of the over- 
running (warm advection) often associated with the XL sector of this synoptic scale 
event. 
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APPENDIX E. MAUI, HAWAII - JULY 1988 



A. METEOROLOGICAL DATA 

Appendix E includes profiles from the July 1 9SS Kihei, Maui rawinsonde dataset. 
A Yaisala DigiCORA MW 1 1 rawinsonde system was used as the data acquisition sys- 
tem. Launch site for this session was located at 20.75° latitude, 156.43° longitude. The 
initial 2 km were sampled every two seconds (S-10 m). Above 2 km, the sampling rate 
was once every ten seconds (40-45 m). Flight number and launch times are provided in 
Table II. 



Table 11. JULY 1988 KIHEI. MAUI RAWINSONDE LAUNCHES 



FLIGHT 

NO. 


TIME (L'TC) 


DATE (DAY MON YR) 


1 


0703 


7 July 1 9SS 


7 


1537 


7 July 19S8 


3 I 0536 


S July 19SS 


4 


1612 


8 July 1988 


5 


0605 


9 July 19SS 


6 


152~ 


9 July 19SS 


7 


05.56 


10 July 1988 


S 


1532 


10 July 1988 


9 


0550 


11 July 1988 


10 


0926 


11 July 1688 



July 00 - Flight 1 




Fig. 62. Kiliei. Maui Rauinsondt* Data: 0703 U1C 7 July 1988 
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July 88 — Flight 1 

M ixing Ratio (g /kg) 




fig. 63. Kiltei. Maui Rauinsoiuic Data: 0703 UTC 7 July 1988 
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July 88 - Flight 1 



Wind Direction (deg) 




Fig. 64. Kiliei, Maui Raninsonde Data: 0703 UTC 7 July 1988 
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July 80 - Flight 2 




Fig. 65. Kihei, Maui Rawinsonde Data: 1537 UTC 7 July 1988 
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July 08 - Flight 2 
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fig. 66. Kiliei, Maui Rnuinsonde Data: 1537 UTC 7 July 1988 
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July 08 - Flight 2 

Wind Direction (deg) 




Fig. 67. Kihei. Maui Rauinsoiide Data: 1537 UTC 7 July 1988 
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Jul} 7 80 



Flight 3 




Fig. 68. Kiliei, Maui Rawinsonde Data: 0536 UTC 8 July 1988 
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July 80 - Flight 3 
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Fig. 69. Kihei. Maui Rauinsonde Data: 0536 UTC 8 July 1988 
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July 00 - Flight 3 
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Fig. 70. Kiliei, Maui Ran insonde Data: 0536 UTC 8 July 1988 
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July BO - Flight 4 




Tig. 71. Kilit'i. Maui Rauiusoude Data: 1612 UTC 8 July 1988 
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July 08 - Flight 4 
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Fig. 72. Kihei, Maui Ran insonde Data: 1612 UTC 8 July 1988 
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July 88 - Flight 4 
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Fig. 73. Kiliei. Maui Rnninsonde Data: 1612 UTC 8 July 1988 
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July 00 - Flight 5 




Fig. 74. Kihei, Maui Rawinsonde Data: 0605 UTC 9 July 1988 
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July 88 - Flight 5 

Mixing Ratio (q/kq) 

0 2 4 6 0 10 12 14 16 




fig. 75. Kiliei, Maui Ran insonde Data: 0605 IJTC 9 July 1988 
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July 80 - Flight 5 

Wind Direction (deg) 




Fig. 76. Kiliei, Maui Ran insonde Data: 0605 UTC 9 July 1988 
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July 08 - Flight 6 




77. 



Kihci, Maui Rauinsonde Data: 1527 UTC 9 Ju)v 1988 
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Fig. 78. Kiliei, Maui Ra»insonde Data: 1527 UTC 9 July 1988 
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July BO - Flight 6 
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rig. 79. Kiliei, Maui Rawinsonde Data: 1 527 UIC 9 July 1988 
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July 00 - Flight 7 




Fig. 80. Kiliei, Maui Raw insonde Data: 0556 UTC 10 July 1988 



129 



July 08 - Flight 7 
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Fig. 81. Kiliei, Maui Rau insonde Data: 0556 U1C 10 July 1988 
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July 88 - Flight 7 

Wind Direction (deg) 

0 50 100 150 200 250 




Fig. 82. Kiltei, Maui Rawinsonde Data: 0556 UTC 10 July 1988 
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July 08 - Flight 0 




rig. 83. Kihci. Maui Rauinsonde Data: 1532 UTC 10 July 1988 
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July 88 - Flight 8 
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Fig. 84. Kihei, Maui Rauinsonde Data: 1532 UTC 10 Julj - 1988 
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July 80 - Flight 8 

Wind Direction (deg) 

0 50 100 150 200 250 




Fig. 85. Kihei. Maui Rauiusoude Data: 1532 UTC 10 July 1988 
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July 00 - Flight 9 




86. Kiliei, Maui Rauinsonde Data: 0550 UTC 11 July 1988 
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Fig. 87. Kiliei, Maui Rauinsonde Data: 0550 UTC 11 July 1988 
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July 88 - Flight 9 
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Fig. 88. Kilici. Maui Rauinsonde Data: 0550 UTC 1 1 July 1988 
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July 00 - Flight 10 




Fig. 89. Kilici, Maui Ran insonde Data: 0926 UTC 11 July 1988 
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July 88 - Flight 10 
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Fig. 90. Kiliei, Maui Rauinsomle Data: 0926 UTC 11 July 1988 
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July 00 - Flight 10 
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91. Kiliei, Maui Ran insonde Data: 0926 U7C II July 1988 
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